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Abstract. 

The scalar meson /o(1370) - indicated in particular in the low energy pp — > 3 body reactions - is a crucial element in 
certain schemes of the scalar meson spectroscopy including glueballs. The most definitive results can be obtained from elastic 
and inelastic nit phase shift analyses using the constraints from unitarity where the discrete ambiguities can be identified 
and resolved. We reconsider the phase shift analyses for k + k~ — > k + k~ , % % , KK, r\r). While a clear resonance signal 
for /o(1500) in the resp. Argand diagrams is seen in all channels above a large "background" from /o(600) there is no clear 
signal of a second resonance "/o(1370)" in this mass range in any reaction, at the level of ~10% branching ratio into %n. 
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ROLE OF / (1370) IN SCALAR MESON SPECTROSCOPY 

There is a particular interest in scalar meson spectroscopy since the lightest gluonic meson is expected within QCD 
with these quantum numbers. The QCD predictions for the glueball mass and a possible mixing with qq states is still 
controversial but the mass falls into a range of 1000 - 1700 MeV. Therefore it is important to identify the low mass 
states with their spin parity, say below 2 GeV, and to classify them into flavour multiplets and glueballs according to 
their production and decay properties. Different scenarios have been proposed in the past for these lowest multiplets. 

According to a popular scheme the states / (600)/cr, K*(900)/k, /n(980), ao(980) are grouped into the lightest 
nonet either built from qq or qqqq. The next higher multiplet of qq type would include ao(1450) and #*(1430). 
For completion of the nonet two isoscalar states are needed but three are found in this mass range according to the 
Particle Data Group (PDG) [1]: / (1370), /o(1500) and f Q ( 17 10). This offers the possibility to introduce the scalar 
glueball which together with the two isoscalar qq states would mix into the observed three scalar /qS. There are various 
realisations of such a mixing scheme first considered in [2]. 

A critical input in this scheme is the very existence of the three isoscalars, so that there is one supernumerous 
state for a usual qq nonet. However, the state /o(1370) cannot be considered as firmly established. This state appears 
with different parameters in the fits and the PDG quotes mass and width in the ranges M = 1200 — 1500 MeV and 
r = 300 — 500 MeV. There are various modes quoted as "seen" in the listing, but not a single number is considered 
to be established among the Branching Ratios or ratios thereof despite many analyses over the years. Quoted results 
from different experiments are often in conflict. Moreover, we do not see a convincing direct signal, neither a peak 
above a low background nor a clear resonance effect in an energy independent (point by point) phase shift analysis. 
The experimental status of /o(1370) is very different from the one of the nearby /o(1500) where the analyses have 
converged and now there are well defined mass and width parameters with errors <10 MeV and five well established 
branching ratios. Also there are clear direct signals from resonance bands in the Dalitz plots of 3 body decays and also 
from phase shift analyses - what is lacking for /o(1370). 

There may be another problem with the above classification scheme in that the mass of the glueball of a full QCD 
calculation is lower than the often quoted mass around 1600 MeV which refers to QCD lattice calculations in quenched 
approximation. A recent full (unquenched) calculation in fact suggests a lower number, around 1000 MeV [3], which 
is below the above 3 mass states, but more work is needed for definitive conclusions. Likewise, in the QCD sum rule 
approach the lightest scalar glueball is found around a mass of 1000 MeV [4]. 

In view of the problems with /o(1370) a classification scheme without this state has been proposed some time ago 
[5]. It has been argued that the activity in this mass region be related to the high mass part of the broad state called 
now /o(600) /a. Then, the states / (980), a (980), a (1450), A"* (1430) are considered to form the lightest qq nonet 
while the broad low mass state is considered to be of a gluonic nature. The "Breit Wigner mass" of the broad state 
where the phase shift passes through 90° is found around 1000 MeV, to remember we also write /o(600 — 1000); the 
width is of size similar to the mass. 



The existence of / (1370) therefore plays a crucial role for the interpretation of the low mass scalar spectrum with 
or without glueball. We look for an answer from a detailed phase shift analysis of elastic and inelastic nn scattering. 

PHASE SHIFT ANALYSIS OF %% SCATTERING 

A resonance in a scattering process corresponds to a pole of the scattering amplitude as function of the complex energy 
E or s = E 2 . Such a pole 1 / (so — s) can be visualized as a closed circle in the complex amplitude plane for physical 
real s; in case of a moving background a circular structure remains although with some distortion. Examples for the 
movement of the isoscalar elastic nn scattering amplitude in the complex plane ("Argand diagram") obtained from 
global fits to various nn data with circles from the various resonances including /o(1370) can be seen, for example, in 
Refs. [6]. 

In the presence of sizable background the existence of a resonance can be demonstrated if the complex partial wave 
amplitude T( is reconstructed from data in a sequence of mass bins and gives evidence for a resonance circle. With 
sufficient statistics this procedure can establish a resonance even without explicitly fitting parametric amplitudes. 



Data on nn — > nn scattering 

Such data have been obtained from the analysis of reactions like np — > nnn by isolating the one-pion-exchange 
contribution. The measurement of the angular distributions and their moments allows for a decomposition of the nn 
amplitude into partial waves up to some descrete ambiguities classified by the "Barrelet-zeros". The unknown overall 
phase can be fixed by assuming a Breit-Wigner form of the leading resonances p(770), /2(1270) and p3(1690). 
Phase shift analyses for n + n~ — > n + n~ became available from the CERN Munich experiment [8] through different 
analyses: 1. CERN Munich-I [9] with a global K matrix fit and subsequent bin-by-bin analysis with one favoured 
solution; 2. CERN Munich-II [10] with a fully energy independent analysis resulting in 4 different solutions, two 
found acceptable. This analysis included all correlations between the angular moments and is carried out in larger 
mass bins resulting in smaller errors for the phase shifts in comparison to CM-I; 3. Similar results with 4 solutions in 
an analysis of CM-I data have been obtained by Estabrooks and Martin [11]. The two acceptable solutions by CM-II 
are shown in the "Argand-diagrams" of Fig. 1. Also shown is the 7 = 2 component used in the present analysis which 
becomes inelastic and moves inside the circle above 1300 MeV (see also Ref. [12]). 

An extended experiment with polarised target by the CERN-Krakow-Munich experiment [13] has allowed decisive 
tests of the underlying exchange models. A phase shift analysis has found a unique solution for nn mass M < 1700 
MeV but with large errors for the parameter 77$ forM > 1400 MeV. Remarkably, the results for the phases 8® from all 
the above analyses agree for M < 1400 MeV. For our present study we take the phase shift results from CM-II which 
have the smallest errors. 

There are also results on the reaction n + n~ — > n°n° obtained from the corresponding np reaction [14, 15]. The 
GAMS data [14] have sufficiently small errors above 1400 MeV to be compared with n + n~ data. The reconstructed 
n°n° amplitude as well as the 7 = component are shown in the rightmost panel of Fig. 1 where the 7 = amplitude 
is found largely inside the unitarity circle as it should be [12]. The results on the 7 = nn amplitude obtained from 
charged and neutral pion reactions should agree. One can see that this is true indeed for the solution ( - + -) (leftmost 
in Fig. 1) at least at the qualitative level: both solutions describe a circular motion of the amplitude with peak near 
1300 MeV and with an additional small circle evolving towards the center near 1500 MeV. The two amplitudes are 
not exactly in agreement but slightly displaced which we relate to systematic uncertainties coming from the 7 = 2 
amplitude badly known in the inelastic region and from the precise form of the line shape of the leading resonances 
used to determine the unknown overall phase. 

The second, small circle is not seen in Sol. ( ). It represents the state /o(1500) with a diameter of about 0.3 

which corresponds to the elastic width and is in good agreement with the PDG result (34.9 ±2.3)%. This agreement 
is remarkable as our number comes from the imaginary part of the elastic nn amplitude whereas the PDG value is 
derived from the measurement of all channel cross sections. 

On the other hand, there is no sign of any additional resonance circle in the 1300 MeV range neither in the CERN 
Munich nor in GAMS data which would reflect the existence of /o(1370) in the way suggested from the global fits 
in [6]. Assuming we should have seen a circle of size 1/3 of the one from /o(1500) we can give an upper limit 
B(/o(1370) ->• nn) < 10%. 




FIGURE 1. Argand diagrams (Im r/ vs. Re Tf) for isospin waves Sj: The two solutions for So ( — I — ) and So ( ) for 

7l + 7C~ -> n + K~ (CM-II [10]) and for n + x~ -> 7r°7r° (GAMS [14]) and its component So (from left to right). Also shown on the 
left panel is the isotensor component S2. Numbers for Tin masses in MeV. 



Data on kk — > KK scattering 

Again such inelastic processes can be obtained from the corresponding np reactions. The experiments with highest 
statistics are: 1. Argonne experiment on n~p — > K~K + n [16]; with primary energy of 6 GeV and KK mass range 
M K g < 1 .55 GeV. Besides the amplitudes with K, A\ , A2 exchanges contributing to KK production there are now also 
(0, B, Z exchanges and the S wave can be produced with 1=1. Then the complete analysis yields 8 different phase 
shift solutions; they are all eliminated but one. 2. Brookhaven experiment on n~p — > K°K°n [17]; this experiment 
with 3-4 times the statistics of other K°K° experiments and with a 23 GeV pion beam has the extended mass range 
M K g < 2.4 GeV. There are even partial waves only and the analysis yields just two solutions below 1.7 GeV related by 
a sign ambiguity. For these reasons we compare here with the results from this experiment. Both experiments [16, 17] 
agree roughly on the S wave magnitude but indicate a difference in phase at low mass M K g < 1 .2 GeV. 



5-MATRIX FIT TO ISOSCALAR S- WAVE 



Next we fit the discussed data from CM-II and BNL with a parametrisation of the coupled channel S matrix in the range 
1 .0 < M nj[ < 1 .7 GeV together with data from CM-I from the range 0.7 < M %% < 1 .0 GeV including the resonances 

a//o(600- 1000), /o(980), / (1500), [/o(1710)] (1) 

where we also take into account the tail of /o(1710) with parameters from PDG but free coupling strength. 

A popular ansatz is based on the addition of several K matrix poles, but the relation to the T matrix poles is not 
obvious and one cannot easily switch on or off a particular resonance. Also it is not obvious why K matrix poles should 
be added, in the elastic case K = tg8. We prefer using a model where the phases are added in the elastic region which 
is obtained if for two resonances the corresponding S matrices are multiplied S = S1S2 which may be generalised to 
one inelastic and one elastic resonance (or background) according to [18]. This leads for the combined T matrix to 
the well known expression T = T2 + T\ exp(2;§2). We generalise this to an Ansatz for inelastic resonances: the broad 
"background" from /o(600 - 1000) interfering with narrow resonances / (980) and / (1500) MeV 

5 = S1S2S3 • • • (2) 



We assume here simple Breit-Wigner forms for the resonances ;' coupling to channels k, i 



Sl = l + 2m where T iM = ^ _ - _ (3) 



We perform a common fit of this parametrization to the experimental data in nn and KK channels for the parameters 
inelasticity 77 and phase 8 for the amplitude of %% -> ^^by T\\ — (77exp(2;'<5) - l)/2i and to |5| 2 and <j) for %% KK 




FIGURE 2. Common fit to TtK and KK data: results for Argand plot (Im Sq vs. Re So), phase shifts 8® in degrees and inelasticities 
r/Q for S-wave k + tz~ 



by T\2 — \S\exp(i<p). We include in the fit the 3 (+1) resonances in (1) each one with parameters M, T and decay 
fractions x, = Ti/T, T = £T, where appropriate line shapes for the channels %%, KK and An are included. The An 
channel is kept to satisfy unitarity but without experimental input. 

The results for elastic nn scattering are presented in Fig. 2. It can be observed that the main features of the data are 
reproduced: the slowly moving background whith phase shifts passing 90° near 1000 MeV, interrupted by the strong 
variation from /o(980) and a second rapid variation near 1500 MeV. One can see from the Argand diagram that the 
fitted circle near 1500 MeV is a bit rotated away from the circle in the data which leads to a shift of the dip in rj away 
from the experimental value. We assume that this is a shortcoming of our rather simple parametrisation with minimal 
resonances and no further "background". However, it does not appear justified to us to include an additional resonance 
to describe the effect which would generate an additional circle. The amplitude (both tj and 8) behaves smooth near 
1300 MeV without any sign of an additional /o(1370). 

Moving on to the KK channel in Fig. 3 we note again a good overall agreement with a rapid variation near 1500 
MeV in both magnitude and phase of the amplitude and a smooth behaviour otherwise. The discrepancy in magnitude 
near threshold cannot be resolved without disturbing the agreement with magnitude (t]q) in the elastic channel and the 
fit represents a compromise. The inelastic amplitude has been fitted including an arbitrary overall phase. 




FIGURE 3. Common fit to TtK and KK data: results for Argand plot, phase shift and magnitude |S| 2 for S-wave tc + tc — > KK. 

There are also data on the reaction n + n~ — > 7777 [19] for which the Argand diagram has been constructed [5]. Again, 
there is a nice "small" circle near 1500 MeV above a slowly moving background but no indication of a second circle. 
These data are not yet included in the fit. 

The shortcomings of the above fits in certain details suggest some improvements. First, it should be noted that 
the parametrisation (2) is to be considered as a phenomenological ansatz; it has correct limits into elastic or mixed 
elastic-inelastic channels but in its general form it is unitary but not necessarily symmetric as it should be. While our 



application is similar to the mixed case a formula with all constraints fulfilled would be desirable. It may be interesting 
to compare with a K matrix ansatz as well despite its other shortcomings. Also one may test resonance formulae 
beyond the minimal form (3) and the role of possible additional background. 

SUMMARY AND DESIRABLE FUTURE STUDIES 

Using data from n + K~ elastic and charge exchange scattering and some information on the %% 1 — 2 channel we are 
able to find a unique solution for the isoscalar %% S wave. The respective Argand diagram clearly shows the small 
circle from /b(1500) consistent with known elasticity above the broad background from /o(600 — 1000) (its directly 
observable Breit-Wigner mass is ^1000 MeV). A simple parametric fit reproduces the main features of the data in 
terms of the broad /o(600 — 1000) and the narrow /o(980) and /o(1500) while there is no direct signal from /o(1370) 
which yields the limit B(/ (1370) nn) < 10%. 

The main difference to the earlier studies with positive evidence for /o(1370) is our technique of using an energy 
independent analysis of the scattering data. In this way the extraction of partial wave amplitudes and their parametric 
representation can be discussed separately. Note that the phase shift data as in Fig. 1,2 together with those from the 
other partial waves come from a very good fit to the original experimental moments with small % 2 < 1 in all bins. 
Therefore there is very little room for additional resonance circles. A global fit to the original data (moments, Dalitz 
plot) typically has a rather large x 2 /ND > 2 and this may be due to a simplified parametrisation (as in our fit) or also 
to an inadequate amplitude structure, the latter case is avoided in performing a point by point phase shift analysis first. 

It will therefore be desirable to perform energy independent analyses with other data of high statistics. Such 
developments have started for 3 -body final states recently [20], noteworthy the recent amplitude determination in 
D s — > 3n by BaBaR [21]. Such data may be precise enough to resolve in a model independent way small circles, i.e. 
resonances with small couplings, such as the suspected fo ( 1 370) . Similar studies could be performed for pp — > 3 body, 
which needed in addition an assumption on nucleon spin couplings. While such analyses cannot determine absolute 
branching ratios from a particular channel as we can do in our analysis, they can help to establish the existence of such 
states like /o(1370) if a positive resonance signal is detected. 
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